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ABSTRACT 


The spherical means of the solutions of a linear partial differential 
equation Lu - f in a conical region are studied. The conical region is 
bounded by a surface generated by curvilinear ? surfaces. The spherical mean 
is the average of u over a constant ^ surface. The conditions on the linear 
differential operator, L, and on the orthogonal coordinates ?, rj, Q are 
established so that the spherical mean of the solution subjected to the 
appropriate boundary and initial conditions can be determined directly as a 
croblem with only one space variable. Conditions are then established so 
that the spherical mean of the solution in one conical region will be pro- 
portional to that of a known solution in another conical region. Applica- 
tions to various problems of mathematical physics and their physical in- 
terpretations are presented. 
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1 . 


Introduction 


In a study of the solutions of the diffraction of an incident wave by a 
cone of solid angle C5 with zero normal derivative thereon it was shown [1] 
that the resultant value at the vertex of the cone is equal to 4rT/(4rT-rj) 

i 

times the value of the incident wave at the same point in absence of the cone. 
This relationship was obtained by rederiving Kirchhoff's formula [2] for the 
wave equation from Green* s formula, for the vertex of a cone. Corresponding 
results can be obtained for the Poisson equation and the diffusion equation. 
Although the results are different for different differential equations, they 
can be summarized in a unified statement; 

"The value at the vertex of the cone is equal to that of an extended 
three-dimensional problem without the cone divided by the ratio of the local 
solid angles of the domains i.e., by the ratio (4tt-0) /(4tt) . In the extended 
problem the inhomogeneous term and the initial data are assigned to be zero 
in the interior of the cone. In case that the normal derivative on the 
surface of the cone is not zero, the prescribed data should be redistributed 
and added to the inhomogeneous term for the extended problem." 

This unified statement suggests that it may be proved directly from 
certain properties common to the differential operators of those equations 
without making use of the Green f s theorem and Green *s function for each 
differential equation individually. Such a general approach can be used 
when the value at the vertex is identified as the limit of the mean value 
over the local spherical cap exterior to the cone. Therefore, the original 
plan of presenting the special analyses for each of these differential 
equations separately is abundoned in favor of the present general approach. 
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In the next section conditions on the linear differential operator L and 
the curvilinear orthogonal coordinates §, r|, C are established so that the 
"spherical" mean of the solution of Lu - f in a domain D subjected to 

i 

appropriate boundary and initial conditions can be determined directly as 

i 

the solution of a problem with only one space variable. The domain D is 
bounded by a "conical” surface S c generated by the ? -coordinate lines and by 
two truncating constant ^-surfaces and S 2 . The "spherical" mean u(§,t) is 
defined as the mean value of the solution u over a constant § -surface in the 
closed region D. Conditions are then established so that the "spherical" 
mean of the solution in D will be proportional to that of a known solution in 
a different region D. For the special case that the surface at 5 = 5 ^ 
degenerates to a point, the vertex of the "cone", the relationship for the 
spherical means is also valid for the value u at the vertex. These state- 
ments are the essence of the Theorems I and II and the two corollaries in 
Section _2. The aforementioned unified statement on the value at the vertex 
of a cone is included as a special case of Corrollary 2 . 

The direct determination of the "spherical" mean u(£,t) provides some 
important information about the solution and also provides a test of the 
accuracy of the numerical approximation to u itself (see, e.g. [ 3 ], [ 4 ].) 

Sections 3 to 6 present the application of the theorems to various 
physical problems. They can be read independently. 

Section ^ applies the theorems to biharmonic equations governing the 
deflection of a circular sector plate subjected to a stationary or moving 
load. The mean deflections and in particular the deflection at the tip are 
identified as the symmetric deflections of a circular plate with equivalent 
loads . 
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In section 4, the theorems are applied to unsteady three-dimensional 
diffusion problems in a wedge and in a domain bounded by two cone 3 
with a common vertex. The demonstrations of how to obtain results directly 
from the theorems due to the particular geometries of the domains are, of 
coirse, not restricted only to the diffusion problems. The third example 
shows their applications to the diffusion of pollutants in a moving medium. 

In section the theorems are applied to the diffraction of an incident wave 
by a cone of solid angle 0. The spherical integral of the resultant wave is 
equal to that of the incident wave without the cone. The value at the vertex 
is equal to the value of the incident wave intensified by the factor 4tt/(4tt-0) 
as stated in Theorem III. If the incident wave is reflected and diffracted 
at least by part of the surface of the cone for all time, the primary wave 
which is known say for t < 0 is composed of the incident wave and those 
reflected and diffracted waves. The spherical integral of the resultant wave 
for t > 0 can again be related directly to that of the primary wave as stated 
in Theorem TV. Their applications are demonstrated by three examples. They 
also illustrate how to apply the theorem to bodies which are composed of cones 
of finite length by making use of the finite speed of propagation. 

In section ^6, the theorems are applied to the reduced wave equation. 

In particular, for the diffraction of a plane wave or a cylindrical wave 
by a cone of infinite length, the Sommerfeld radiation conditions are not 
applicable. These theorems provide a necessary condition for the resultant 
wave . 
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2 . 


Direct Determinability of the Mean-value Over a Curvilinear Coordinate 
Surface 

Let u(x,y,z,t) be the solution of a linear partial differential equation 

(2.1) Lu = f(x,y,z,t) in D for t > 0 

where x,y, and z are the Cartesian coordinates and t denotes the time variable 
Let g, r], C he a set of curvilinear orthogonal coordinates. The element of 
arc length is given by 

(2.2) dx 2 + dy 2 + dz 2 = a 2 dg 2 + 0 2 d r} 2 + y 2 d£ 2 

The domain D is bounded by two constant g -surfaces and and a cylindrical 

surface S ^ generated by g -coordinate lines with f 2 > g > (see Fig. 1). In 

the g, n, Q space, the domain D is a right cylinder parallel to the g-axis 
with base area B in the r)-£ plane. The initial and boundary conditions will 
depend on the operator L. It will be assumed henceforth that these conditions 
and Eq. (2.1) form a well posed problem. To be more specific, Eq. (2.1) is 
assumed to be second order in g and the standard boundary conditions of the 
third type [5] will be imposed, i.e., 

(2.3) a jU + (-1)^ b^ u^ = hj (Tl,C,t) on for j = 1,2, 

where a and b. are non-negative functions of t and a. + b. is greater than 
J J J J 

zero. The required boundary condition on S c will be specified later in 
Theorem I, The initial conditions at t - 0 are 
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(2.4) 



(x,y,z) 


in D for k = 0,1... K-l 


where K is the order of the differential operator L with respect to t. It 

will be assumed, throughout this paper, that all the given delta f , h and q 

j n k 

are piece-wise continuous in their domains of definition, that a, 0, Y process 
continuous first derivatives and that the solution u has continuous first 
derivative in the closed region D i.e., D and its boundary. 

The generalizations of the initial conditions (2.4) and of the boundary 
conditions (2.3) on when Eq. (2.1) is of any order in ? will be discussed 


after Theorem I. 

The "spherical" mean v(f,t) of a function v(5 jTlj £,t) is defined as the 
mean value of v over the constant 5 surface in 5, i.e.. 


(2.5) 


V(!,t> - J/ V dS/ JJ dS 

S s 

= JJ V aydndC / JJ SYdridS 

B B 


Since the "spherical" mean involves only one space variable, it would be 

of interest to know under what conditions on the operator L and the curvilinear 

coordinates and with what type of boundary conditions on the surface S will it 

c 

be possible to determine the "spherical" mean u(?,t) directly prior to the 
solution u(§,ri,£,t) itself. - 

The linear differential operator L and the curvilinear coordinates are 
assumed to fulfill the following four conditions: 
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( 2 . 6 ) 


(i) L = L x + L 2 
where involves only g and t. 

(ii) 1*2 u can be written as the divergence of a vector A in a 
constant g surfaces, i.e., 


(2.7) 


L ’ “ - div A = ^ [ 2- (ccyA^) + <08A r )] 


where A and A- are linear combinations of u and its partial derivatives and are 

T) Q 


(iii) The operator commutes with the averaging operator ( ) 


identified as the T]- and ^-components of the vector A with A^ = 0. 

(iii) The op 
defined by Eq. (2.5), i.e., 

( 2 . 8 ) 


(Ljii) = LjU . 


When involves differentiation with respect to g, the dependence of the area 
element of the surface S on g should be separable, i.e. 


(2.9) 0 Y = V(S) M(T1,0 • 

(iv) The distance between two constant ^-surfaces is a 

constant, i.e., 


(2.io) a - a(g). 

On the surface S , the required boundary condition is to prescribe the 
normal component of A i.e. 
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( 2 . 11 ) 


A n = g(?,a»t) on S c 


where a is the arc length along the intersection of S, and S . If there are 

S c 

more than one boundary condition on in the well-posed problem Eq. (2.11) 
has to be a linear conbination of them. 

The divergence theorem is applied to the vector function A for the 
volume AD in D bounded by two constant % surfaces S_ and S, + Af. By 

9 3 

utilizing Eq. (2.1), conditions i) and ii) and the boundary condition of 
Eq. (2.11), the result is 

A? JJ aSy a x u - f) dr|dC = i? / a g (§,a)d{r 

S ? * S ? 

Due to condition (iv) , a(£) can be moved outside the integral. After dividing 
the equation by the volume AD, the equation becomes 


( 2 . 12 ) 


where 

(2.13) 


II 

ft 


a iU ) BydridC 


If SydndC 

S ? 


g (§,t) 

3S. 




g (I .o’, c) da / JJ Bv<ir|d£ 


'5 


The leftside of Eq. (2.12) is the ’'spherical" mean of u which is equal to 
u ° n account of condition (iii). Equation (2.12) then becomes the 
differential equation for u. The boundary conditions and the initial conditions 
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for u are supplied by averaging Eqs. (2.3) and (2.4). These results are 
summarized in the following theorem. 

Theorem I 

When the differential operator L and the orthogonal curvilinear coordinates 
fulfill conditions (i) to (iv), the "spherical” mean u(?,t) of the solution of 
Lu = f in D subjected to the boundary conditions of Eq. (2.3) and (2.11) and 
the initial condition of Eq. (2.4) is governed by the following equations: 

(2.14) D.E. L 1 u(?,t) - f(§,t) + g(e,t),for < 5 < ? 2 and t > 0, 

(2.15) B.C. a^u + b (-1) ^ = h^(t)^at % ~ for j = 1,2, 


(2.16) I.C. -2_u = (?) ^ at t . 0 for Jj, < ? < ? 2 and k = 0,1, ...K-l. 

d 1 


From physical point of view, the operations leading to Eq. (2.14) is 
equivalent to setting up the "conservation" equation for a thin shell in 5 
between surfaces and with thickness aA?* If f is interpreted as the 

intensity of a spatial source distribution in D and g as a surface distribu- 
tion of sources over S^, the inhomogeneous term in Eq. (2.14) for the 
"spherical " mean represents a spatial source distribution independent of r\ 
and Q with intensity equal to the "spherical" mean of the spatial distribu- 
tion f plus a redistribution of the source on S c with constant § uniformly 
over the surface S_. 
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V 


It should be noted that in Theorem I, the operators on u in the initial 
conditions and in the boundary conditions at § = and ^ are tlie same as 
those for u because these operators are independent of r; and Q. It is evident 
that if the operators on u in Eq. (2.4) are generalized to be any linear 
operator independent of rj and Q f then the same operators will appear in the 
initial conditions on u. Similarly, the boundary conditions (2.3) on and 
can be generalized with total number of conditions equal to the order of 
Eq. (2.1) in provided that those boundary operators on u are linear and 
independent of T) and Q. The same operators will again be applied to u in its 
boundary conditions on ? - ^ and in Theorem I. 

The geometry of the surface S c or that of the base area B in the r|-£ 
plane appears implicitly in the definition of the inhomogeneous terms of 
Eqs. (2.14) and (2.16) and has no effect on how to construct the solution u of 
those equations. The solution u will not be changed if there is an inter- 
change of the surface source distribution g with the spatial distribution f 
so long as the sum f + g; remains unchanged. In particular, all the surface 
source g can be replaced by an increament of f in D by the amount *g with u 
remaining the same. In physical problems, the inhomogeneous term f can contain 
point sources in D, therefore, the surface source distribution can be absorbed 
in f directly. In this respect, the discussions can be simplified by imposing 
a homogeneous boundary condition on S , i.e. 

(2.17) A = 0 

n 

* 

There are cases that u (§,Tj,C,t) which is the solution of the same set of 
equations for u but in a different domain D , is given or can be obtained 
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—if ic 

easily. The "spherical 1 ' mean u(C,t) can be computed directly from u by the 

definition of Eq. (2.5). It is therefore of interest to relate the "spherical" 
* 

mean in D to that in D. 

Domain D is also bounded by two constant 5 surfaces and S 2 with 
end §2 respectively and by a cylindrical surface generated by ^-coordinate 
lines with > 5 > with base area B* in the rj-£ plane. 

From condition iii) on the curvilinear coordinates, the ratio of the area 
of constant - \ surface in D* to that in D becomes 

(2 * 18) . Jj^ydridC / JJ Qydr\d Q = JJ |_l(n,QdndC / JJ M(T1,Odr]dC = A 

B B* 


where A is a constant. 

It will be assumed that the supports in inhomogeneous terms in the 

differential equation and in the initial conditions lie in D fiD*, and the 

supports of the inhomogeneous terms in the boundary conditions on S. lie in 

B f) B*. The solution u* in the domain D* is then defined as the associated 

solution of u in D if they fulfill the same differential equation, the same 

boundary conditions on and S 2 , and the same initial conditions, and the 

homogeneous condition A ~ 0 on S for u and A * 0 on S * for u* where A* is 

n c n c — 

the vector ^A with u replaced by u*. The last condition implies either 

g = g* = 0 or that the surface source distribution g has been absorbed in f. 

-* * 

Let v designate the "spherical" mean of v over a constant 5 surface in 

* --- 

D . The mean values of the inhomogeneous terms in the differential equations 
initial conditions and boundary conditions on S . for u will be A times the 
values for u, e.g. , 
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(2.19) 


f * = A f 


By theorem I, u* obeys the same equations as u with the exception that all 
the inhomogeneous terms change by a factor A, therefore, 

(2.20) u* (5,t) = A u(5,t) or u(?,t) = u*(?,t)/A . 

This result is stated as follows: 

Theorem II 

The "spherical” mean u(?,t) of the solution u in D is equal to "spherical" 
mean u*(5,t) of the associated solution u* in D* divided by the ratio of the 
area of the constant 5 surface in D to that in D*. 

By using the definition of Eq. (2.5), Eq. (2.20) becomes 

(2.21) JJ u(§,Tl,C,t)dS = JJ u(5,n,C,t)dS 



This relationship restated in the following Corollary is useful when one of 
the areas, say that of S*, is not finite. 

Corollary I 

The integral of u(§,T},C,t) over a constant 5 -surface in D is equal to 
that of its associated solution u* in D*. 

*L. Admissible Orthogonal Coordinates 

Conditions iii) and iv) i.e. Eqs. (2.9 and 2.10), impose some 
restrictions on the curvilinear orthogonal coordinates. It is shown in an 
NYU Report [6] that Eq.(2.10) requires the ^-coordinate lines to be straight 
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lines then the unit tangent vector 5 is independent of When % is identified 
as the arc length, a(?) equals to 1. It is then shown in [6] that the constant 
5-surfaces and the orthogonal coordinates, which are consistent with Eq. (2,9), 
have to be one of the following three types: 

a) 5 i s independent of T) and £. The constant f surfaces are parallel 
planes say normal to x-axis. The variable 5 can be identified as x while 
T}»C can be identified as y and z or a pair of orthogonal coordinates in the 
y-z plane. The domain D is a right cylinder parallel to x-axis (see Fig. 2a). 

b) 5 ^- s dependent only on T\ and Q. The constant 5 surfaces are coaxial 
circular cylindrical surfaces. If the axis is identified as the z-axis. The 
orthogonal coordinates can be identified as the cylindrical coordinates (p,8,z) 
respectively. The domain is bounded by the cylindrical surfaces p = 5^ and 

p = ?2 an( * t ^ ie "conical" surface S c which can be defined as S(z,9) = 0, 

(see Fig. 2b). 

c) 5 depends on both rj and £. The variable 5 is identified as the 

radial distance r to the center and T] and £ can be identified as the spherical 
angles 9 and cp respectively or as two orthogonal coordinates on the spherical 
surface. The domain D is a cone truncated by spherical surfaces of radii 5^ 
and 52* ( see 2c). 

2 

Under type c, the area S of a spherical cap in D is r f] where ft is 

r c c 

the solid angle of the cone. When = 0 the spherical cap at r = 5^ 
degenerates to the vertex of the cone and the spherical mean value u(0,t) 
becomes the value at the vertex of the cone. Theorem I will therefore relate 
the value at the vertex of the cone D to that at the center of a spherically 
symmetric problem without the conical surface S^, i.e.. 
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( 2 . 22 ) 


u (r = 0, 9,cp,t) = u(0,t) 


The boundary condition(s) at r - ^ = 0 will be the specification of the 
behavior of u and/or its derivatives with respect to r as required by the 
order of Eq. (2.1) with respect to The same condition(s) will be imposed 
on u at r » 0. 

When D* is another conical domain of solid angle ft^ the ratio of the 

area of spherical cap S in D to that of S* in D* is 

r r 





When the required boundary condition for u on S^, Eq. (2.11), is 
homogeneous. Theorem II yields that the spherical mean of the solution 
u(r,0,(E,t) in D is equal to that of u in D divided by A which is the 
ratio of the solid angles, i.e., 

(2.23) u (r,t) = J*(r,t) / (n*/n ) 

c c 

— 

At r - 0, the mean values u, u are the values of u and u* respectively at 
the vertex. The following statement is now evident. 

Corollary 2 

The solution u at the vertex of a cone D with 0 < r < ^ are equal to 
the associated solution u* at the vertex of cone D* divided by the ratio of 
their solid angles. 

For the special case that D* is the whole sphere, r < with ft* = 4tt, 
the associated solution u* is a solution inside the sphere in absence of a 
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conical surface, u* at the center of the sphere is related to the solution 
u and the vertex of the cone D directly, 

(2.2A0 u(r = 0,9,cp,t) = u* (r = O,0,co,t) 

\ c 

\ 

This equation is equivalent to the statement in Section 1 which summarizes 
the results obtained by the applications of Green's formulas and the 
appropriate functions to several equations of mathematical physics. 

It should be pointed out here that conditions for the direct deter- 
mination of the spherical mean in Theorem I are not necessary. For example, 
the operator involving the variables, f| and Q can be removed by means of 
the Stokes Theorem if L^u can be expressed as the ? -component of the curl of 
a vector. Details of the alternate approach can be found in [6]. They will 
not be presented in this paper because in most physical problems, the 
governing equations are usually derived from the conservation of certain 
quantities and the representation of u as the divergence of a surface 
vector is usually expected. 

In the next sub -section it will be shown that several equations of 
mathematical physics fulfill the conditions stated in Theorem I. 

B) Equations of Mathematical Physics 

For many equations of mathematical physical, namely the Poisson equation, 
the unsteady diffusion equations, the wave equation, etc., the differential 
operator involving the space variables is the Laplacian operator, i.e. 

(2.25) L = div grad + 

m 

where involves differentiations with respect to time only. 
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To compare with the conditions on L, the vector A in Eq. (2. 7) can be 
identified as the projection of grad u on the surface, i.e. 


(2 26} a - — ^ rt + ~ ^ C 

(Z,Z6) - e an n + y a£ & 

where £ a nd £. are unit vectors along the T| and ^-coordinate lines respectively. 
From the definition that Lu = L^u + div A, the operator is [7] 


L„ « 


1 

a$Y 


2u 


8y_ 

a 


b 


+ L. 


When the coordinates are admissible in the sense defined in the preceding 
subsection becomes 


(2 - 27) L x ■ *fr a^ Cvc?) $1 1 + L t 

It is clear that involves only 5 and t and the operator L defined by (2.25) 
fulfills the four conditions required for the direct determination of the 
"spherical* 1 mean. The normal component of A on is equal to the normal 
component of grad u. The required boundary condition on Eq. (2.11), 
becomes 


(2.28) 


g <S,CT,t) 


on S 


In the analysis of diffusion through a moving incompressible fluid, the 
equations for the velocity field £_ is uncoupled from the equation for the 
diffusion process. For a given velocity field, the operator L is of the 

m 

form 
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(2.29) 


L s "f" • V “ A 


Lu can be rewritten as u + V . [qu - vu] . To be consistant with the conditions 

for Theorem I, it is necessary that the ^-component of the velocity is 

independent of Tj and Q and that (£ . n)u-(3u/dn) is prescribed on the surface 

S . The second condition is equivalent to the prescription of the flux of u 
c 

through the surface S^. For the special case that the surface S c is a rigid 
wall, n vanishes on S c and the required boundary condition on S c is the 
prescription of ^u/^n i.e. again Eq. (2.28). 

For biharmonic equations, say 


(2.30) L = Lb + L fc 

the vector A in Eq. (2.7) can be identified as 


(2.31) 


* = 31 c F In Au + 8 L S & an 5 ] 


.1 5u. 
'3 2 3TV 

I SHi 


+ L t 7 ** + Y L ? ( ^2 tf> 1 

where L = -1 — (v ^~=) . From the definition of Lu = L 1 u + div A, it can be 

g v(?) as aS 1 

seen that the operator is equal to (L^J + L fc and involves only t and 5 . The 
required boundary condition on S in Theorem I becomes 


(2.32) A n = |^ Au + n . Q, 0Lg (p|~-) + H • £ Y = 8 (?»<r,t) 

It is obvious that the operator L in general will not be invariant with 
respect to a translation of the origin, see for example Eq. (2.29). However, 
many of the operators do have this invariant property, e.g., Eq. (2.25) and (2.30), 
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and then Theorem I can be used to yield the mean-value theorem for any interior 
point and to develope solutions for initial value problems. 

C) Connections with the Mean-Value Theorem and Initial Value Problems 

In this subsection it is assumed that the operator L fulfills the conditions 
in Theorem I and possesses the property of invariance with respect to a transla- 
tion of the origin. Any interior point P. of D can be chosen as the origin 
and Theorem I can be applied to determine directly the spherical mean u inside 
a sphere lying in D and centered at P^. The spherical mean and in particular 
the value of u at P q . is related to the boundary data on by solving a problem 
with only one space variable. Thus a mean value theorem for Eq. (2.1) is 
established. Of course the boundary conditions on the conical surface will 
not appear in this consideration. 

For an initial value problem, the domain D is the whole space. Any point 
P Q can again be chosen as the origin and Theorem I will relate the solution at 
P Q to the initial data by solving directly the equation for the spherical mean. 
Applications of spherical means to initial value problems in this respect can 
be found in [5] and [8]. 

The theorems in this section and the examples in the following sections 
deal with initial boundary -value problems. Additional examples can be found 
in [6]. 
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3. 


The B1 harmonic Equation 

The biharmonic equation appears frequently in elasticity problems 
namely in the plane stress or plane strain problems [9] and in the deflec- 
tion of thin plates [10]. It is also the governing equation for a slow 
viscous flow [11]. In this section u(r,0,t) will represent the deflection 
of a thin circular sector plate in polar coordinates r,0. The governing 
equation is 


(3.1) 


Mu + 


3 2 u 

3T 7 


= f 


in D • 


The domain D is the circular sector of radius R and angle oC, i.e., o < r< R 
and -°/2 < 0 < a /2. For the two dimensional problem, eq. (2.31) reduces to 


(3.2) 


A = 0 { I 

— — l r 


f- Au + i-1 r iL- ( 1, 
90 3r l 9r 1 r 


— )1 I 

90 


° A 0 


Under theorem I, the two boundary conditions on S c , i.e., on the 
two straight edges 6 = ± a /2, should be sufficient to specify A 0 . As a 
mathematical problem, A 0 can be specified as one boundary condition and 
the mean solution will be independent of the other boundary condition. 

The two boundary conditions to be imposed for the plate will 
appear naturally in the formulation of the governing equations from the 
variational principal [10]. The following terms appear in the line inte- 
gral along the boundary 

-M n + (Q n - 9 M ns 

9n 3s 

where M n , M ns and Q n are the bending moment, torsional moment and shearing 
force and n and s are arc lengths normal to and along the boundary respec- 
tively, The line integral vanishes when one of the following four pairs 
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Of boundary data are prescribed: (i) u and 3u/3n, (11) u and M p , (ill) M n 
and (Q n - 3M ns /3S), (iv) 3u/3n and (Q n - 3M nS /3S). When the boundary condi- 
tions are homogeneous, the first three pairs are known as built-in, simple 
supported and free edge conditions- The fourth pair will be called rotation- 
constrained edge conditions, e.g., the edge is reinforced by a bar with 
torsional rigidity much larger than its flexural rigidity to constrain the 
deflection of the plate. 

Along the straight edges of the sector plate,6 =± a /2,only the 
fourth pair of the edge conditions will yield the data on A g . They are 

(3.3) . Uq = 0 and 

(3.4) 1 (Au) + (1 - v) 3 2 (1 3u ) = o 

r 39 3r 2 r 30 

for 0 < r < R where v denote Poisson ratio in this section. From eq. (3.3), 
eq. (3.4) reduces to 3(Au)/3G = 0. From the definition of A Q in eq. (3.2) 
it is clear that the combination of egs. (3.3) and (3.4) yields 

(3.5) A 0 = 0 

on 0 = ± a /2. The two boundary conditions at the tip of the plate, r = 0, are 

(3.6) u is finite and 

(3.7) lim r 3_ f _J_ __3 , r _3u } 1 = 0 # 

r + o 3r L r 3r k 3r J 

Equation (3.7) implies there is no concentrated load at the tip, otherwise 
the limit will be non zero. The boundary conditions along the circular 
arc r = R can be, 

(3.8) u = 0 and 3u/3r = 0 for a built-in edge , 
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(3.9) u = 0 and 3 2 u/3r + (v/r) 3u/3r = 0 for a simply supported 

edge for |0j < a /2. The boundary conditions are at r = 0 and r = R do not 

i 

involve 0, therefore. Theorem I can be applied to the sector plate to deter- 
mine the mean u (r»t) directly and in particular the value at the tip 
ju(O,0,t) which is IT (0,t). The mean tr in this two dimensional problem 

l 

iis defined as 

m 

(3.10) tr (r,t) =J u (r,S, t) de/a 

r - a/ 2 

The following two examples will deal with the static and unsteady problems 
respectively. 

Example 1. Deflection of a sector plate due to a stationary concentrated 
load . For a concentrated load P applied at a point Cr 0 , 0 Q >in D the inhomo- 
genous term f in eq. (3.1) becomes 

(3.11) f = (P/I) 6(r-r 0 ) 6 (0 - 0 o )/r 

where I is the flexural rigidity of the plate. With a built-in edge along 
the circular arc (eq. 3.8) and rotation-constrained edges along the straight 
edge (eqs. 3.3, 3.4), the mean deflection u (r) is governed by the equation 


(3.12) 


(r i- r ,r = (P/ 1 ) <5(r-r 0 ) /(m) 

^ dr dr 


The boundary conditions at r = 0 and r = R are given by eqs. (3.6), (3.7) 

and (3.8) with u replaced 17 . The solution is 

(3.13) u (r) = 4 ^{( r2 + r 0 2 ) ^g ^o + (R + r 2 ) (1 - )/2 

R R 

+[r 0 2 - r 2 ) + (r 0 2 + r 2 ) log (r/r 0 )] H(r - r 0 )J 

* 

where H is the Heaviside unit step function. The deflection at the tip 
of the sector plate is 
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(3.14) u(0,e) = u (0) = _P_ [r 0 2 log »o. + (R 2 - r 0 2 ) /2] 

4<xl IT 

The mean deflection is the same as the symmetric deflection [10] 
of a circular plate with built-in edge at r = R and a uniformly distributed 

i 

lioad along a concentric circle of radius r 0 with total strength P(27r/a). 

| 

Example 2. Deflection of a sector plate due to a moving load 

For a concentrated mvoing load of strength P, the inhomo- 
geneous term f (r,t) is again given by eq. (3.11), however, r 0 (t) and 
0 O (t) are now given functions of t describing the motion of the load. 

If the boundary conditions on the sector plate are the same as those 
in example 1, the mean deflection tt (r,t) is governed by the equation 

(3.15) r |_) 2 y + v « £• <S [»" - r Q (t)l/ (ra) 

The boundary conditions atr = 0 and r = R are eqs. (3.6), (3.7) and 
(3.8) with u replaced by u. If the initial condition for u is u = u t = 0 
at t = 0, then the same conditions hold for IT, i.e. 

(3.16) u (r,0) - u t (r,0) = 0 

It should be noted that eq. (3.15) does not involve 0 Q (t). Hence, the 
mean deflection TT (r,t) is independent of the circumferential movement of 
the load and remains the same if the load P moving along any radial line 
with instantaneous position r Q (t). XT (r,t) is equal to the symmetric 
deflection of a circular plate with built-in edge at r = R, with same 
initial conditions and subjected to a load of total strength P (27r/a) 
uniformly distributed along the circle r = r 0 (t). The solution for 
the symmetric deflection can be constructed by making use of the eigen- 
functions for the symmetric modes given in[12]. Procedures for the con- 
struction of the symnetric unsteady solution are described in [13] and [14]. 



4. The Unsteady Diffusion Equation 

If uCx,y,z,t) represents the temperature variation in 
a stationary medium or the mass fraction of species in a diffusion 

problem, u obeys the simple diffusion equation 

(4.1) Lu = Au - C -2 u^_ = f(x,y,z,t) in D for t > 0. 

where C 2 the diffusivity constant. In addition to the boundary 
conditions of eqs.(2.3) and (2.11), there is one initial condi- 
tion , , 

(4.2) u(x,y,z,0) = q 0 (x,y,z) in D 

Two examples are presented to illustrate how to select 
the curvilinear coordinates for a given domain and how to apply 
the Theorems and corollaries. These suggestions, which rely on 
the geometry, of the domain, are also applicable to other problems 
of mathematical physics. Of course, the final explicit formulas 
will be different for different governing equations. The first 
example illustrates the use of Theorem II to obtain explicit 
solutions along the edge of a wedge, and those for the planar, 
the cylindrical and the spherical integrals in the wedge. The 
second example shows the use of Theorem II for the construction 
of the spherical mean in a domain bounded by two cones with a 
common vertex. The third example shows how to apply Theorem I 
to study the diffusion of pollutants in a moving medium and in 

particular to obtain directly the total flux of pollutants going 
upstream. 

Example 1. Three dimensional heat transfer problem in a wedge 
Figure 3 shows that the domain D is a wedge with angle 
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y. The edge is chosen as the z-axis. One face of the edge is the 
x— z plane (<$> = 0) » Due to the method of superposition and the 
Duhamel*s principal, it suffices to consider the simple case 
of homogeneous equation and boundary data while the initial data 
is a point source located at a point P(x ,y ,0) with strength Q, 
i.e.,f=0,g=0 and 

(4.3) q (x,y,z) = Q 6(x-x , y-y ,z) 

Q 11 

with x J = p cos <f> , y = p sin cj> , p >0 and 0 < <f> < y. 

ii iii ii i 

In absence of the wedge, the solution u* of the homo- 
geneous heat conduction equation in the whole space under the 
initial condition of eq.(4.3) is [15] 

(4.4) u*Cx,y,z,t) = Q(4irC 2 t)" 3/2 exp{ -r 2 /(4C 2 t) } 

where r is the distance from (x,y,z) to the location of the 

l 

source . 

In term of cylindrical coordinates p,<f>, and z, z can 
be identified as the £ coordinate and the planar integral of the 
solution u in the wedges at constant z is equal to that of u* 
in the whole space by Corollary 1, i.e. 

co y 

I (z,t) = / pdp / d<J> ,z ,t) 

Poo 

oo ao 

= / dx/ dy u*(x,y,z,t) 

_oO 

(4.5) = Q( 4ttC 2 1 ) ~ 1 / 2 exp[ -z 2 /(4C 2 t) ] 

If now P is -identified as £ . A constant £ -surface is 
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£ cylindrical surface with radius p. Corollary 1 gives 
immediately the cylindrical surface integral of u, i.e., [ 7] 

00 |j 

I (p,t) = / dz / Pd<J> u( p ,<£ ,z ,t ) 

| c _00 0 

1 _ 2 71 

1 . 

s / dz / Pd4> u*(p cos p sin <t>,z,t; 

_co 0 

= Q p(2 C 2 t) -1 I 0 [PP i /(2C 2 t)] exp{ -( P 2 

< 4 - 6 > + p 2 )/4C 2 t} 

Since any point A(0,0,z o ) on the edge can be identified 
as the vertex of a cone with solid angle 2 n the point A can be 
chosen as the' origin of the spherical coordinates (rj.co). The 
value of u at A is therefore related to u* by corollary 2, 
u(0,0 > z o ,t) = — u* ( 0 ,0 ,z 0 ,t ) 


(4.7) = (2tt/u) Q(4irC 2 t)" 3 /2 exp{-[p 2 + ( z ^ -Z q ) 2 ] / ( 4C 2 t ) } 

The spherical integral of u with A as the center is also related 
to that of u* as follows 


I (r ,z ,t) - v 1 f sin 9d0 / d<fu - r 2 f sin 6d6 / d* u* 
s 0 o o » o 


or 

I (r,z c ,t) = Q(4TtC 2 t) -I / 2 ( r /r p ){exp[-(n - >\) 2 ] 

s 

(4.8) - exp[-Cn + n i > 2 ]} 


where ^l 2 = r /(4C t) 


D 2 = r 2 /(4C 2 t) 
i P 
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is the distance from the 


and r p = [p* + (z^)]^ 

point source P to A. 

i 

1 Explicit formulas are obtained for the solution along 

i 

the edge in eq.(4.7) and for the integrals of the solution on 
a constant z-plane, and a cylindrical surface with z-axis as 
its axis and with radius p and on a spherical surface of radius 
r centered at a point (0,0,z Q ) on the edge in eqs . (4 . 5 ) , (4 . 6 ) 
and (4.8) respectively. The only missing information about the 
solution is its dependence on the variable (J> . The success of 
obtaining so many explicit results from the corollaries is due 
to the special geometry of the domain, a wedge, and is, therefore, 
not restricted to diffusion problem only. 

Example 2 . Temperature field in a solid bounded by two insulated 
conical surfaces with a common vertex. 

As shown in fig. 4, the inner and the outer conical 
surfaces are nonintersecting. Their solid angles are Q and 
respectively with > fi . With the boundary condition of 
3u/3n = 0 on the two conical surfaces, theorem II and the 
corollaries remain applicable with the solid angle of the domain 
D equal to - Q . Based on the reason stated in the first 

O J 

example, it suffices to consider only the simple case of a point 
source of strength Q at a point P in D at the instant t = 0 . 

The solution u* in the whole space without the conical 
boundaries is again given by eq.(4.4). The temperature u at 
the common vertex A is again given by the corollary 2 . 
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u(A,t) = 


4 it 

- 0, 


u*(A,t) 


_ 2 
JLl. 


(4.9) 


= ( 


4tt 


- 4C 2 t 


«o ~ a , } (4TTC 2 t) 3 / 2 


where r = | PA| is the distance from the point source to the 

common vertex A. The integral of the temperature over the 

spherical surface in D with radius r and centered at A 

is again given by the expression on the right side of eq.(4.8). 

When the conical surfaces are of finite length £ 

the domain D is bounded in addition by an insulated spherical 

surface S with radius £ . The associated solution u* should be 
2 2 

replaced by the solution in an insulated sphere of radius £ 

2 

due to a point source at P at the instant t = 0. The value at the 

vertex u(A,t) and the spherical mean can be related by Theorem I 

to a spherically symmetric problem given in [isJ • 

The axially symmetric temperature field between two 

co-axial circular cones is investigated by Oberhettinger and 
Dressier [ 16] Series solution is constructed for the isothermal 
boundary condition; u = O . For case of adiabatic boundary 
condition, 3u/9n = 0, only a brief outline of the method is 

presented. The explicit results in this example would be use- 

- 4 

ful to check the numerical results of their analysis. 
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Example 3. Diffusion of pollutants in a moving medium . 

To be more specific, the conical domain D will be an open 
channel connecting two large reservoirs. The channel is assumed to 
be of constant depth, h. The domain D is defined by 

(4.10) < l < 5 a »-ct/2 < n < a/2 and -h < ? « z <0 

For a constant area channel £ * x and n = y; a is the width of the 
channel. For a divergence channel, £ = r, and n = ©; a is the diver- 
gent angle of the channel. 

The flow is assumed to be incompressible and the velocity field 
can be represented as 

A 

q = I V ( 5 ,/d 

is the velocity at the section £ ■ A = 0 for a constant area 
and A * 1 for a divergence channel. 

The governing equation for the diffusion of a pollutant in the 
is [15]. 

la. + vr 5 i) X M -c 2 au = c 2 f . 

at Vl 5 ' 35 

The boundary condition on S c i.e., on the wetted walls of the 
and on the free surface is 

9u = 0 , 

an 

The boundary conditions at the upstream and downstream ends of 
the channel will be 

u = 0 at £ * Ci and £ = £ 2 


(4.11) 

where V 
channel 

channel 

(4.12) 

channel 

(4.13) 
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because of the large sizes of reservoirs. The inhomogeneous term f 
represents the source of pollutant in the channel. For unsteady prob- 
lems an initial condition, say u = 0 at t - 0, shall be added. 

j If the pollutant is dumped at a point £ Q , n 0 » Z Q in D at a 

rate Q(t), then f can be written as 

(4.14) f = Q 6(5 -5 0 ) 6(n - n 0 > «(z - *<,) IK 

The mean solution u (£, t) obeys the equation 

(4.15) || + V( |/ | - c *l0 + 5' X H-J = C=Q «(? -<*)/(**«) 

The boundary conditions are 

(4.16) u = 0 atC = and £ 2 

The mean solution u (£,t) can be obtained directly and is inde- 

pedent of n 0 and z Q of the point source. From u (£,t), is obtained 

inmediately and is the quantity of interest in this problem. In particular, 

u at £i» and £2 represent the rate of the pollutant entering the upstream 
C 

and downstream reservoirs. 

Explicit results will now be given for the steady problem to show 
the difference between a constant area channel and a divergent channel. 

The steady solution of eqs. (4.15, 4.16) for a constant area 
channel (X = 0) is 
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SV'V (,k -«‘ (£ ° ' 5,)k ) [1 - « ■ E °V»" “ 3 

j*V for €„<£<£, 

-(C 2 - e)k -(C - €i)k / [1-e -JtX J 

% [1-e J tl-e ° 1/ 

av 


for Z 2 <Z<Z q 

where k = V/C 2 and X = £ 2 - £j. The percentage of the pollutant reaching 
the upstream reservoir is 


(4.17) 


aC 2 


U^Si) 

Q 


*v e 


e "(«o * 5 i )k 
-(? 0 - Ci)k 


tl-e 

as k oo 


-kX 

e 0 ) k i/ u-«» i 


The percentage decreases exponentially in the distance from the source of 
the pollution to the upstream station of the channel. 

For a divergent channel {X = 1, V > 0), the steady solution is 


Re Re Re 

M) "w-'Sndt^ -'J 


Re 


-f - *J H f5-«o>J 


where Re = V^i/C 2 . The percentage of the pollutant going to the upstream 
reservoir is 


Re 


Re 


Re 


tr'fCi) ci a „ (W U u - i . ( £i» r i- (g o / . eO . . 1 

(4.19) Q (ej/g,) «e_i ^i*, 1 t 1- U,/??) Re J 

As Re -► « f the percentage decreases as the Re-th power of the ratio §j/£ . 


When V < 0, the flow is reversed from £ - £ 2 to £ * £i through 
a convergent channels Equations (4.18) and (4.19) remain valid for J/ <0, 
however, eq. (4.19) defines now the percentage of pollutant going to the 
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downstream reservoirs and 

-Re 

(4.20) u * JL- ft. ^ i . (.Si) as - Re -+ °° 

Q $2 

The percentage going to the upstream reservoir can be computed from 
$2 a. l or can be seen from (4.20). It again decreases as the 
(-Re)-th power of the ratio (£ 0 /£ 2 ) as -Re -* °° in contrast with the 
exponential decay law in a constant area channel. 
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5 . The Wave Equation 


The propagation of acoustic waves or that of a single 

component of electro or magnetic waves is governed by the simple 

\ 

wavej equation 


(5,1) Lu = Au - C~ z u = f (x ,y , z ,t ) 

tt 


where C is the speed of propagation. The usual two initial 
conditions are the prescription of u and u^_ at a given instant 
. say t = 0 . 

In diffraction problems, it is usually done to 
define an incident wave u^x,y,z,t) and to designate t = 0 as 
the instant the incident wave front hits the body. In other- 

(i) 

words, the support of u does not intersect the body for t <_ 0 . 

The solution u(x,y,z,t) for t > 0 D # i.e., in the domain exterior 

to the body is a solution of -the wave equation, eq. (5.1) subjected 


to the initial conditions of 

(5.2) u(x,y,z,0) = uHx,y,z,0) 
and 

(i) 

u (x,y,z,0) = U (x,y,z,0) 

(5.3) t t 


and the boundary condition 


(5.4) 


u = 0 on S 

n c 


where S is the surface of the body. When the body is a cone of 
c 

solid angle ft, its exterior angle ft , which is 4tt - ft, is the 

. c ... Ci> . 

solid angle for the domain D. Since the incident wave u is a 
solution the wave equation in the whole space without the cone 
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for t > 0 with the same initial conditions of eqs.(5.2) and 

g * v 

(5.3), u can be identified as the associate solution u* in the 
whole space D* and corollaries 1 and 2 yield 

<i> 

u (o.o,o,t) 

with S* in D* Equation (5.6) is the result obtained previously 

' - r 

[1] by a different method. 

If the resultant wave is written as the sum of the 
incident wave and a secondary wave, which is the reflected and 
the diffracted waves, i.e., 

<i> 

( 5 . 7 ) u = u + w 


(5J5) ff uds = ff fcis 

! S S * 

1 r r 


and 

(5.6) 


u(o,o,o,t) = 


4tt 

4tt - ft 


eq.(5.5) becomes 


(5.8) 


ff w dS = ff u^s 

s s*-s 

r r r 


These equations can be summarized in the following statement. 


Theorem III 

When an incident wave is reflected and diffracted by a 
cone of solid angle after the instant t = 0, the spherical inte- 
gral of the incident wave cut off by the cone is redistributed 
exterior to the cone as the integral of the reflected and the 

diffracted waves. The redistribution of the incident wave from 

* 

a spherical surface of angle 4 tt to that of 4 tt - Q intensifies 
the value at the vertex by the ratio 4tt/(4tt ~ ft ). 
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If an instant (say t = 0) cannot be found such that 
the incident wave is neither reflected nor diffracted by the 

core 9 theory I^^cannot be applied since the initial data can- 

■«* 

not be stated in terms of the 'incident wave alone. In other- 
words 9 Theorem II cannot be employed to relate the solution u 

in D tf i.e., exterior to a conical body D^to u^ in the whole 

space without the body. However t it still can be used in diffraction 

problems to relate the solution u to a known solution u* in D* which 

£a exterior to another conical body I>£ with the same vertex as stated 

ih the following theorem: 

Theorem IV 

If the support of the solution u* outside a conical body 
Djj does not intersect the conical body for for t£Q and the ^ur- 

It 

face of the cone contains the part of the surface of where 

the incident wave is reflected or diffracted for t£0, then the 
solution u in D , fulfills the following relationships 
(5.9) u 5 u for t i 0, 


u dS = // 

S* 

r 


and 

(S.il) = < 4* V? u (0 ’ t} » for arL t 

where S and S* are the spherical surfaces outside the conical bodies 
r r 

D and D* respectively „ The spherical surfaces are centered 
c c 

at point 0, the vertex^ and the solid angles of the conical 

bodies are Q and Q* respectively. 

The two solutions u and u* differ from each other for 
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t>0 because that the reflection and diffraction will then take 
place over different surfaces. However, u and u* have the same 
initial data at t = 0, therefore, Theorem II is applicable and 
eqs.(5.1Q) and (5.11) are valid for all t. 

Theorems III and IV are special adaptations of Theorem II 
for diffraction problems. Their applications are demonstrated in 
the following three examples. The first example shows a direct 
application of Theorem IV to a trihedral and the explicit formulas 
become useful for tie more complicated problems in the next two 
examples. The second and the third examples deal with bodies 
which are composed of cones of finite length. These examples 
illustrate how to apply the Theorems to each conical component 
for a finite interval due to the finite speed of propagation. 

Example 1. Diffraction by a trihedral. 

Figure 5 shows an incident wave which is diffracted 

all the time by one edge (0A) of a trihedral with solid angle 

ft. The primary wave u ^bill pass over the vertex or one of the 

other two edges OB and 0E at the instant t = 0 . The given 

primary wave u^ p ^represents the Incident and the diffraction waves 

advancing over a dihedral of angle X formed by the two plane OAE 

and OAB . The solution u exterior to the trihedral will differ 
(p) 

from u only for t > 0. 

If the planes OAE and OAB are extended beyond the edges 
OB and OE to form a wedge in place of the trihedral, u P will be 
the solution also for t > 0. Point 0 which lies on the edge of 
the wedge can be considered as the vertex of a conical body of solid 

angle 2y. With u* identified as and ft* as 2y, 'Theorem IV yields 
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(5.12) 


u( 0 ,t ) 


4tt - 2y 
4 tt - ft 


(P) 

u (0 ,t) 


and 


“IT 


2 TT -X 


(5.13) // u dS = r 2 f sin 9 d0 / d<£ u (p)(r , 0 , t) 

S o o 

r 

where r is the distance to the vertex 0 , S r is the spherical sur- 
face outside the trihedral, 9 = 0 is the edge OA and cp = 0 is the 
plane OAB. 

Example 2 . Diffraction of a plane pulse by a rectangular block 
(a two dimensional problem) 

As shows in Figure 6 # a plane pulse of strength e hits 
the first corner 0 of the block at the instant t = 0 at an angle 
X with respect to the horizontal side OA. The pulse is reflected 
by the two sides of the corner and the influence of the vertex is 
confined by the circle of radius Ct . The solution outside the 
circle and behind the reflected wave is 2e. The solution inside 
the circle is a conical solution [171 in two variables, u c (x,y), 
with x = x/(Ct) and y = y/(Ct) . Before the construction of the 
conical solution, the "spherical" integral and the value at the 
vertex 0 will be related directly to the incident wave . 

For the two dimensional problem, cylindrical coordinates, 
p,0 with the vertex 0 as the origin will be used. The “spherical " 
integral will become circular integral of radius p. Theorem III, 
that is, eqs. ( 5.5 and 5 . 6 ) yields the circular integral of the 
solution u due to the diffraction by the first corner at 0 

3 IT/ 2 2 TT 

/ u (p, 0 ,t) pdQ = / u^(p, 0 ,t)pd 6 „ 

0 1 0 


(5.14) 


^e(2irp) for p<^Ct 

e[2.7T-2 cos -1 (Ct/p)]P for p > Ct > 0 
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t > 0 


and the value at the vertex 0 , 

, 2 tte 4 

(5.15) u (0,t) = — 7 = - e 

l 37T/2 3 

Equation (5.14) is valid only when the circle of radius p does 
not receive the diffracted waves from the adjacent corners. 

This condition creates an upper time limit for the validity of 
eq . ( 5 .14) . It is 

(5.16) Ct < L (1 + sin X ) - p 

For equation (5.15), the time limit is given by eq.(5.16) with 
p = 0. For the convenience of explanation, the width L of the 
block is assumed to be shorter than its height H and the incident 
pulse arrives at the corner A(L,0) before at the corner B(0,H) 
i.e,, L sin X < H cos X is also assumed. 

If the conical solution u c (x,y) is constructed by the 
method of C17J, it defines the solution u-^ for the domain exterior 

to the first corner at 0. Theorem ivcan be used to relate the 
diffracted wave at the adjacent corner A, to u^ when t > L sin X/C. 

For a finite duration, corner A can be considered as the vertex 

of a wedge. In terms of the cylindrical coordinates p and 0 , 

i l 

the domain D exterior to corner A is -it/2 <0 < it and p >0. 

i l 

The associated domain will be the space above the line AO, and 

its extension through A as shown by the dotted line in Fig . 65 

i.e. 0 < 0 < - tt and p >'0. For the associated domain D*, there 

1 i 

is no diffraction at the corner A, therefore, the associated 

solution can be identified as u for t > L sin X/C. The circular 

^ * 

integral of the solution u^ exterior to the corner A is given 

by that of u , 

1 
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TT 1? 

/ u (p ,0 ,t )d9 - / u (p ,0,t)d0 

— TT / 2 211 1 0 1 1 1 

for 2L + L sin X - p >_ Ct > L sin X. The ratio of the angles 
of the domain D* and D at the vertex A is tt/(3tt/ 2) = 2/3. The 

for L >_ Ct > L sin X 

for 2L + L sin X >_ Ct L 

These equations can also be obtained by applying the results of 
example 1. It is clear that if the solution u^ after the 
diffraction at corner B is constructed by the method in C 18 3 
Theorem IV will be useful for the subsequent diffraction of 
at an adjacent corner. 

Example 3. Diffraction of a plane pulse by a cube. 

Figure 7 shows a plane pulse of strength e hitting the 
first corner 0 of a cube at the instant t = 0. The edges QA, 

OB and 0E are chosen as the x -, x and x - axes respectively. 

12 3 

The directional cosines of the normal to the plane pulse are 
n , n and n . They are all non-negative and the smallest one 

12 3 

is n j . Before the arrival of the plane pulse at any other corner, 

i.e., t < n L/C, the solution u (x.,t) will remain unchanged if 
l i i 

the cube is replaced by the three dimensional corner at 0 with solid 
angle Q = 4tt/ 8 = ir/2. Equations (5.5) and (5.6) yield the 
spherical integral of u outside the corner, 

tt it . 

// u (x.,t)ds = / d0 / u 1 r 2 sin 0 d4> • 

S 1 1 o -tt 

r 

er 2 4 < 7T tor r <_ Ct 

erLr + Ct] 2 tt for r > Ct > 0 



value at the vertex A is 


2 4 

T(2e) = -e 


f 3 “ 

U ( A,t ) = \ 3 

2 Uu ( l 5 

Q c Ct 5 


0) 
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and for Ct < (1 + n )L ~ r and "the value at "the vertex 0 
, Hire 8 

u^COjt) = 4 _ ir ^ 2 = j e for (1 + n )L > Ct > 0 

The upper time limits are defined by the arrival of 
the diffracted waves from an adjacent corner. Along the edges 
and outside the domain of influence Q f the vertex i.e., r > Ct , 
the value is same as that along the edge of a 90° wedge, i.e., 

u (x . ,0 ,0 ,t ) = (4/3)e for Ct/n. > x. > Ct, 

ii l i 

i = 1,2,3. 

To obtain the complete solution u (x,y,z,t) it is 

necessary to apply the method of £173 for conical solutions of 

two variables outside the characteristic sphere r = Ct and the 

method of [4] for the three dimension conical solutions inside 

the sphere . After the determination of the solution u for 

i 

the diffraction by the corner 0, Theorem ^ and in particular, 
the results in example 1 can be applied to the solution of 
diffraction of u^ by the adjacent corners. 

When t > n^L/C, the wave front of u advances over the 
edges AB r and AE T and the vertex A. A can be considered as the 
vertex of a trihedral with solid angle ft = it/2, and 0A as the 
edge of a dihedral with angle x = tt/2 . The results of example 1 
i.e., eqs . (5.12) and (5.13) are applicable until the arrival of 
diffracted waves from other corners. The result V 'are 

(5.17) U 2 <A,t) = 6/7 u^Ajt) . 

and 
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(5.18) 


ff (u - u )dS = 0 

S , 2 1 

r 

where u is the solution after the diffraction of u by the 
2 1 

comer at A, S_^i is the spherical surface outside the trihedral 

with radius r f to its vertex A and the extension of u outside 

i 

its domain gf definition is equal to zero. Similarly, eqs. (5.12) 
and (5.13) can be applied for any point on the edge AB f (or AE ! ) 
before the arrival of the diffracted waves from the vertex A or 

B*. 

The applicability of eqs . (5.12) and (5.13) to vertex 
A to obtain eqs. (5.17) and (5.18) does not depend on the facts 
that u^ is a conical solution with respect to 0 and that 0 and A 
are vertices of a cube. It is obvious that the results of example 
1 can be applied to relate the solution u* of diffraction of a 
wave by a vertex A of a polyhedron to the solution u of the 
subsequent diffraction of u* by an adjacent vertex B of the poly- 
hedron. If the polyhedron is convex, the solution u can be written 
as u* + u \ The extension of u* outside its domain of definition 
equals to zero. There is a discontinuity in u* and also in u^ J but 
their sum is continuous. Theorem IV or Eqs. (5.17 and 18) yield 
CdJ 

(5.19) // u ds = 0 
S 

r 

and 

u <i> (B,t) = u(S,t)-u*(B,t) 

2 x - 8 

(5.20) 4rr jj" u *(B,t) 
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where x is the angle of the dihedral with edge AS and 0 is the 
solid angle of the polyhedron at the vertex B. Of course, eqs . 
(5.19 and 5.20) will remain valid only before the arrival of 

I 

additional diffracted waves from the adjacent corners. 
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6 ■ The Reduced Wave Equation 


For periodic solutions of wave equation v(x,y,z)e , 
the amplitude function v is governed by the reduced wave equation 

(6.1) Av + k 2 v = f in D 

where k 2 = c 2 u) 2 . 

If the domain D is bounded and in particular a truncated 
cone (£ 2 < *) the governing equation for the spherical mean v(r) 
is an ordinary differential 

(6.2) ( v f r 2 ) * + k 2 r 2 v = f 

There is no difficulty to obtain the spherical mean v(r) for 

£ < r < £ and the value at the vertex of the cone when £ = 0 . 

1 — — 2 1 

When the domain D is unbounded, there is the well- 
known difficulty of what is the proper boundary condition at in- 
finity. This difficulty does not appear for the solution of an 

initial boundary value problem of the wave equation. 

For the diffraction of a incident periodic wave with 
amplitude v lJ by a obstacle, the solution v can be represented in 
the form 

«i> , 

(6.3) v = v + w 

where w is the secondary wave, the diffracted and reflected waves. 
When the obstacle is finite, the boundary condition at infinity 
is the well known Sommerfeld radiation conditions t 19 J 

♦ 

(6.4) w = Q(r -1 ) and + i k w = oCr** 1 ) 

The conditions can be deduced by treating the periodic solution 
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as limit of an aperiodic solution for example the principle of 
limiting absorption or limiting amplitude £20 J . 

If the obstacle is unbounded, say a half plane, a wedge 

or a cone, and the incident wave itself fulfills eq.(6.4), 

S <i) . 

e.g.j v is induced by a point source, or a source distribution of 
bounded support, the radiation conditions of eq.(6.4) for the 
secondary wave remain valid [ 21 J . 

If the incident wave itself does not fulfill the Sommer- 
feld radiation conditions, for example, the incident waves are 
plane waves or cylindrical waves, the secondary wave due to 
an unbounded obstacle may not obey the Sommerfeld 
conditions. These difficulties are mentioned in [22] and 
suggestions are made to separate the part of w which fulfill 

the Sommerfeld condition from the second part, which does not. Explicit 
determinations of the second part are demonstrated £227 a plane and 

a half plane and can be done for a wedge . For a conical body with 

arbitrary cross-section, the determination of the second part in w 
is not obvious. 

In the first subsection, the incident wave is assumed 
to fulfill the Sommerfeld conditions and so does the secondary 
wave. The corollaries of Theorem II are applied to relate 
the spherical mean of the resultant solution to that of the 
incident wave and to show that their relationships are the same as 
those of unsteady solutions of the wave equation. These relation- 
ships for the unsteady solution^ will now bg assumed to be valid for 
the periodic solution including the case that the incident wave does 
not fulfill the Sommerfeld conditions and yield the necessary condi- 
tions for the solution at infinity in the second subsection. 

42 


•pT — T 


The incident wave fulfills the Sommerfeld conditions 


Since v^ } and w both fulfill the Sommerfeld conditions, 

so does the resultant solution vin the domain D outside a conical body 

* 

of solid angle ft. Without the body, i.e. in the whole space D , 

• * • 

then there is no secondary wave and the solution v is equal to 

v a> and fulfills the same homogeneous boundary condition on r = £ + 00 9 


l.e. eq. ( 6 .4) # 


. Corollaries I' and II are therefore 


applicable, and the results are: 


( 6 . 5 ) 


// v dS : // v dS 


(6.6) v(0 ,0 ,0 ) = vH 0,0,0) 

A 

where is the spherical surface with radius r from the origin 
located at the vertex of the conical body and S r is the part of 
the surface outside of the body. Equations (6.5) and (6.6) are 
identical to the relationships for eq.(5.5) and (5.6). In terms 
of the secondary waves w, eq.(6.5) becomes 


(6 .„ // ■ {£ s 

r r r 


ci), e 
v dS 


This is identical to eq.(5.8) and has the same physical meaning 
stated in Theorem IV, i.e., the spherical integral of the secondary 
wave outside the body recovers the part of the spherical integral 
of the incident wave cut-off by the body. 


B. The incident wave does not fulfill the Sommerfeld conditions 

, T 

In this case, the proper boundary condition at infinity 
for the solution v or the secondary wave w is unknown, therefore. 
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the corollaries cannot be applied. However, if the limiting 
principle is accepted also for the diffraction by an unbounded 
body then the periodic solution ve^ wt can be considered as 

s , , 

thd> limiting solution (t-^-oo) of the wave equation with an incident 

i 

wave which is periodic behind a wave front advancing towards the 
cone at finite time (say t=Q). Equation (5.8) or Theorem IV 
which holds for all time should be valid also for the limiting 
solution. By this argument, eq. (6.7) is reestablished, 


( 6 * 8 ) 


ff w dS = 


ff v^ds 

s£-s r 


and also 


( 6 . 9 ) 


ff wdS 
S r 


= ff 

S*-S 
r r 


v ds 
r 


They become now the necessary conditions for the secondary 
wave w for any r including r -*■ 00 . 

By looking at the Green’s formula for the reduced 
wave equation [ 2 3 the integral over the spherical surface of 
large radius R is 


-kr 

3 /£ ) 

[ // vdS] 4Trr 

S R r = R 


-kR 

[// v ds] 

S r 4 itR 
R 


The terms inside the first square bracket is recognized as the 
total strength of doublets on and that inside the second one 
is that of sources on S^. Equations (6.8) and (6.9) can there- 
fore be interpreted as the following necessary conditions: 

The total strength of sources induced by the secondary 
wave on the spherical surface of large radius R outside the conical 

body recovers that induced by the incident wave on the part of the 
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spherical surface cut off by the cone. The statement holds also 
for doublets. 
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Admissible Orthogonal Coordinates and the Constant 
£ -Surfaces : 

(a) Parallel planes, 

(b) Co-axial circular cylindrical surfaces and, 

(c) Concentric spherical surfaces 






Fig. 3 Three Dimensional Problem in a Wedge 
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Domain D bounded by two cones with a common vertex 








Fig. 7 


Diffraction by a cube 





